Detailed analysis of the tetraquark potential and flip-flop in SU(3) lattice QCD 
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We perform the detailed study of the tetraquark (4Q) potential V4Q for various QQ-QQ systems 
in SU(3) lattice QCD with /3 = 6.0 and 16 3 x 32 at the quenched level. For about 200 different 
patterns of 4Q systems, T4q is extracted from the 4Q Wilson loop in 300 gauge configurations, with 
the smearing method to enhance the ground-state component. We calculate V4Q for planar, twisted, 
asymmetric, and large-size 4Q configurations, respectively. Here, the calculation for large-size 4Q 
configurations is done by identifying 16 2 x 32 as the spatial size and 16 as the temporal one, and 
the long-distance confinement force is particularly analyzed in terms of the flux-tube picture. When 
QQ and QQ are well separated, 14q is found to be expressed as the sum of the one-gluon-exchange 
Coulomb term and multi-Y type linear term based on the flux-tube picture. When the nearest quark 
and antiquark pair is spatially close, the system is described as a "two-meson" state. We observe a 
flux-tube recombination called as "flip-flop" between the connected 4Q state and the "two-meson" 
state around the level-crossing point. This leads to infrared screening of the long-range color forces 
between (anti)quarks belonging to different mesons, and results in the absence of the color van der 
Waals force between two mesons. 

PACS numbers: 12.38.Gc 12.39.Mk 12.38.Aw 12.39.Pn 



I. INTRODUCTION 

The inter-quark force is one of the elementary quanti- 
ties for the study of the multi-quark system in the quark 
model. Since the first application of lattice QCD simu- 
lations was done for the inter-quark potential between a 
uark and an antiquark using the Wilson loop in 1979 
j| , the study of the inter-quark force has been one of the 
important issues in lattice QCD 2] . Actually, in hadron 
physics, the inter-quark force can be regarded as an ele- 
mentary quantity to connect "the quark world" to "the 
hadron world" , and plays an important role to hadron 
properties. 

In addition to the QQ potential 0, 0, ES H| , recent lat- 
tice QCD studies clarify the three-quark (3Q) potential 
Hillll, which is responsible to the baryon structure. 
In fact, our group recently studied the 3Q potential V3Q 
in detail with lattice QCD, and clarified that it obeys the 
Coulomb plus Y-type linear potential [E IE 13 • How- 
ever, no one knows the inter-quark force from QCD in 
the exotic multi-quark system such as tetraquark mesons 
(QQ-QQ), pentaquark baryons (4Q-Q), dibaryons (6Q) 
and so on. 

In these years, various candidates of multi-quark 
hadrons have been experimentally observed. + (154O) 
E3, IH E2, S— (1862) H3 and 6 C (3099) ll| are 
considered to be pentaquark (5Q) states [IE HE] be- 
cause of their exotic quantum numbers, although some 
high-energy experiments report null results [13 ■ A(3872) 




HESSE! and D s (2317) 
tetraquark (4Q) states [24H 2 
from the consideration of their mass, narrow decay width 
and decay mode. 

These discoveries of multi-quark hadrons are expected 
to reveal new aspects of hadron physics, especially for the 
inter-quark force such as the quark confinement force, the 
color-magnetic interaction and the diquark correlation 
[33| . According to these experimental results, it is desired 
to investigate the inter-quark force in the multi-quark 
system directly based on QCD |H HE [H S M, EE 
l4ll |42| , which would present the proper Hamiltonian for 
the quark- model calculation of multi-quarks 0, 13 C3| • 

As for the 4Q candidates, A (3872) was discovered in 
the process B+ K+ + A (3872) ->■ K + + n-n+J/ip at 
KEK (Belle) in 2003 [lal. and its existence was confirmed 
by Fermilab (CDF[l<|^DO|2j3) and SLAC(BABAR)|H. 
£> s (2317) was also found in B-B reaction at T(45) 
resonance at SLAC(BABAR) [^3 and consecutively at 
KEK(Belle) |l3 in 2003. As the unusual features of 
A (3872), its mass is rather close to the threshold of 
D Q (cu) and D°*(uc), and its decay width is very narrow 
as T < 2.3MeV (90 % C.L.). These facts seem to indicate 
that A (3872) is a 4Q state 1241 1251 or a molecular state 
of D°(cu) and D°*(uc) l2rj j Hjga Jg p rather than an 
excited-state of a cc system |3(j, [3l|, |33 • Also, Z) s (2317) 
cannot be regarded as the simple meson of cs, but is con- 
jectured to be a 4Q state from the similar reasons on the 
mass and the narrow decay width. 
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Also in the light quark sector, the possibility of 4Q 
hadrons has been pointed out. For instance, Jaffe 
proposed in 1977 that the light scalar nonet including 
/o(980) and ao(980) can be interpreted as QQQQ rather 
than QQ 46] . Since then, many studies of the scalar 
nonet have been done in terms of the 4Q state 0, E3 • 

As an analytical guiding model for the multi-quark sys- 
tem, the flux-tube picture |4lLl5fll5lll5l l5ll5l l5ll56]] 
has been investigated for the structure and the reaction of 
hadrons, and is supported by recent lattice QCD stud- 
ies H El m m m laa E3, El E3 in this pic- 
ture, valence quarks are linked by the color flux-tube as a 
quasi-one-dimensional object. The flux-tube has a large 
string tension of about a ~ 0.89 GeV/fm, and there- 
fore its length is to be minimized. For the multi-quark 
system, this picture predicts an interesting phenomenon 
of the "flip-flop", i.e., a recombination of the flux-tube 
configuration so as to minimize the total length of the 
flux-tube in accordance with the change of the quark lo- 
cation |5!|,E(|. This process is important not only for the 
structure of multi-quark systems but also for the reaction 
process of hadrons. 

In this paper, we study the tetraquark (4Q) poten- 
tial, i.e., the interaction between quarks in the 4Q sys- 
tem directly from QCD by using SU(3) lattice QCD at 
the quenched level, and investigate the hypothetical flux- 
tube picture for the multi-quark system and the flip-flop 
in terms of QCD. Here, the lattice QCD Monte Carlo 
simulation is the first-principle calculation of QCD and 
is considered as the only reliable method for nonpertur- 
bative QCD at present. We note that lattice QCD is also 
a useful method to select out the correct picture for non- 
perturbative QCD in the low-energy region through the 
comparison with the lattice results. 

The organization of this paper is as follows. In Sect. II, 
after a brief review on the lattice studies of static quark 
potentials, we present a theoretical form for the 4Q po- 
tential based on the flux-tube picture. In Sect. Ill, we 
present the formalism for the 4Q Wilson loop and the 4Q 
potential. The lattice QCD results are shown in Sect. IV. 
In Sect. VI, we compare the lattice QCD results with the 
theoretical form, and discuss the flux-tube picture and 
the flip-flop. Sect. VI is devoted for summary and con- 
cluding remarks. 



II. THEORETICAL CONSIDERATION FOR 
THE 4Q POTENTIAL 

A. QQ, 3Q and 5Q potentials 

To begin with, we give a theoretical consideration for 
the multi-quark potential. From a lot of lattice QCD 
studies [JHHSHI3, * ne sta ^ c QQ potential is known 
to be well described by 
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where r denotes the distance between the quark and 
the antiquark. The first term is considered to be the 
Coulomb term due to the one-gluon exchange (OGE) and 
Aqq denotes the Coulomb coefficient. The second term 
is the linear confinement term with the string tension, 
CTqq ~ 0.89 GeV/fm. 

From the recent detailed studies in lattice QCD with 
(/3=5.7, 12 3 x 24), (0=5.8, 16 3 x 32V(/3=6.0 i 16 3 x 32), 
and 09=6.2, 24 4 ) |E II II II IH IS S E3, the 3Q 
potential is clarified to be the sum of the OGE Coulomb 
term and the Y-type linear confinement term as 
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Here, L m { n denotes the minimal value of the total flux- 
tube length, which corresponds to the Y-shaped flux-tube 
linking the three quarks. In fact, the lattice data of the 
3Q potential, can be well reproduced with only 

three parameters, A^q, 03Q and C3Q. 

To demonstrate the validity of the Y-Ansatz, we show 
in Fig.l the lattice QCD data of the 3Q confinement 
potential V^Q ni , i.e., the 3Q potential subtracted by its 
Coulomb and constant parts, 
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plotted against Y-shaped flux-tube length L m i n , at 
0=5.8, 6.0 and 6.2 in the lattice unit. For each 0, clear 
linear correspondence is found between 3Q confinement 
potential V r 3 c q nf and L m i n as V3Q — <J3QL m i n , which in- 
dicates the Y-Ansatz for the 3Q potential [38L 13^ . liol ] . 

This lattice QCD result strongly supports the flux- 
tube picture for baryons, and the Y-type flux-tube for- 
mation is actually observed in lattice QCD through the 
direct measurement of the action density in the presence 
of spatially-fixed three quarks [33, |33, |32, |4(J, H3 ■ The 
Y-Ansatz for the 3Q system is also supported by recent 
further lattice QCD studies |58l liSSj and analytical stud- 
ies |6fll6l) . 

As for the relation between Vqq and V3Q, we have 
found the OGE result A3Q ~ ^Aqq and the universality 
of the string tension ctqq — er go, which also supp orts the 
flux-tube picture [H Hi EJ M IEIJ5IJ5I Eg and the 
strong-coupling expansion scheme~ f50ll5l| . 

Very recently, the 50 potential is also studied in lattice 

QCD [2lElllEiED,ElE3- lt is wel1 described 
by the OGE Coulomb plus multi-Y type linear potential 
|ill0SE!Ei3,E2- With the minimal length L min of 
the flux-tube, the 5Q potential can be well described as 
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FIG. 2: The tetraquark (QQ-QQ) state. The QQ (QQ) clus- 
ter belongs to 3 (3) representation of the color SU(3). 



For each case, we consider below the theoretical form of 
the tetraquark potential V4Q. 



FIG. 1: The lattice QCD result for the 3Q confinement po- 
tential V 3 C Q nf , i.e., the 3Q potential subtracted by its Coulomb 
and constant parts, plotted against Y-shaped flux-tube length 
imin at /3=5.8, 6.0 and 6.2 in the lattice unit. The clear linear 
correspondence between 3Q confinement potential V 3 c q nf and 
£min indicates the Y-Ansatz for the 3Q potential. 



with (A 5 q, (75q) fixed to be (A 3 q, 03Q) following the OGE 
result and the universality of the string tension. This 
lattice result also supports the flux-tube picture. 



B. Theoretical Ansatze for the 4Q potential 



FIG. 3: The connected tetraquark system. All quarks and 
antiquarks are connected with the single flux-tube, which is 
double- Y-shaped . 



Now, we investigate the theoretical form of the 
tetraquark (4Q) potential V4Q with respect to the flux- 
tube picture, which seems workable for QQ mesons and 
baryons. For the argument of the low- lying 4Q states, 
we consider the 4Q state of ((QQ)g-(QQ)3)i as shown in 
Fig. 2. Here, (QQ)g denotes that two quarks form the 3 
representation of the color SU(3). The meaning of (QQJ3 
is similar. By combining (QQ)3 with (QQ)3, the 4Q sys- 
tem can be constructed as a color-singlet state. We note 
that another possible 4Q system of ((QQ)6-(QQ)g)i i s 
expected to be a highly-excited state, since the attrac- 
tive (repulsive) force acts between quarks, when the QQ 
cluster belongs to 3 (6) representation in a perturbative 
sense, which leads to the 3 diquark picture [33|. 

In the flux-tube picture, the flux-tube is formed so 
as to minimize the total flux-tube length of the system 
for the low-lying state. For the low-lying 4Q system, 
there are two candidates for the flux-tube configuration 
according to the 4Q location. One is the connected flux- 
tube system where all quarks and antiquarks are con- 
nected with the single flux-tube as shown in Fig. 3. The 
other is the disconnected flux-tubes corresponding to a 
"two-meson" state as shown in Fig. 4. Note that the 4Q 
state of ((QQ)g-(QQ)3)i generally includes such a "two- 
meson" state of (QQ)i(QQ)i as shown in Appendix A. 



Qi Q 3 
• o 



• 9 
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FIG. 4: The disconnected tetraquark system, which corre- 
sponds to a "two-meson" state. 



1. OGE plus multi-Y Ansatz for the connected 4Q system. 

For the connected 4Q system, we propose the "OGE 
plus multi-Y Ansatz" as a theoretical form of V4Q from 
the viewpoint of the flux-tube picture. This type of the 
flux-tube configuration is plausible when the distance be- 
tween QQ and QQ is enough long compared with the 
size of these clusters. For such a system, all quarks and 
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antiquarks are linked by the connected double- Y-shaped 
flux-tube as shown in Fig. 3, and the 4Q potential V4Q 
is described by the OGE Coulomb plus multi-Y linear 
potential V C 4Q as 
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with rjj = |rj — r j | and L m ; n being the minimal value of 
the total flux-tube length. Here, denotes the location 
of ith (anti)quark in Fig. 3. 

The first term describes the Coulomb term due to 
the OGE process. Note that there appear two kinds 
of Coulomb coefficients (Aiq, j^Uq) in the 4Q sys- 
tem, while only one Coulomb coefficient, ^4qq or A^q, 
appears in the QQ or the 3Q system. In this defini- 
tion, the Coulomb coefficient A 4 q is expected to satisfy 
A4Q ~ 5^4qq as the OGE results. The brief derivation 
of the OGE Coulomb terms is shown in Appendix A. 

The second term is the linear confinement term with 
the string tension (T4Q, which is expected to satisfy 
04Q ~ Cqq ~0.89 GeV/fm as the universality of the 
string tension. Similar to the 3Q and the 5Q systems, 
the Y-type junction appears in this case, and L m j n is ex- 
pressed by the length of the double- Y-shaped flux-tube 
as shown in Fig. 3. 

In the extreme case, e.g., m = r%s — r 2 A, the 

lowest connected 4Q system takes an X-shaped flux-tube, 
although the energy of such system is larger than that of 
the two-meson state in most cases. 



2. The two-meson Ansatz for the disconnected 4Q system 

For the disconnected 4Q system corresponding to the 
"two-meson" state as shown in Fig. 4, we adopt the "two- 
meson Ansatz" for V4Q. Such a flux-tube configuration 
is plausible when the nearest quark and antiquark pair 
is spatially close and the system can be regarded as the 
"two-meson state" rather than an inseparable 4Q state. 
For such a system as shown in Fig. 4, the total potential 
V4Q for the 4Q system would be approximated to be the 
sum of two QQ potentials in Eq.Q as 



V ±Q = VqQ(ri 3 ) + VQQ(r 24 ) 

= - A Qq(^ + ^) + fT QQ( r i3+ r 24) + 2C Q Q 
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C. The 4Q potential form and the flip-flop 

For the lowest 4Q system, the 4Q potential Vaq would 
be expressed with lower energy of the connected 4Q sys- 
tem or the two-meson system, 



V iQ = min(y c4 Q,^ 2 QQ)- 



(7) 



As a physical consequence of Eq. JJJ based on the flux- 
tube picture, there can occur a physical transition be- 
tween the connected 4Q state and the two-meson state 
according to the change of the 4Q location. This phe- 
nomenon occurs through the recombination of the flux- 
tube and is called as the "flip-flop" . (A popular usage of 
the "flip-flop" may be for the simple flux-tube recombi- 
nation between two- meson states. We here use this term 
as the general flux-tube recombination.) 

The flip-flop is important for the properties of 4Q 
states especially for the decay process into two mesons. 
In addition, the flux-tube recombination between two- 
meson states can be realized through the two successive 
flip-flops between the two-meson state and the connected 
4Q state. Therefore, this type of the flip-flop is important 
also for the reaction mechanism between two mesons. 



III. FORMALISM FOR THE 4Q POTENTIAL IN 
LATTICE QCD 

In this section, we present the formalism to extract the 
static 4Q potential. Similar to the derivation of the QQ 
potential from the Wilson loop, the 4Q potential V4Q can 
be derived from the gauge-invariant 4Q Wilson loop Waq 
as shown in Fig. 5 j3l ifl li^ . 




FIG. 5: The tetraquark (4Q) Wilson loop for the calcula- 
tion of the 4Q potential V4Q. The contours Mi (-1=1,2) are 
line-like, and Li,Ri (i=l,2) are staple-like. The 4Q gauge- 
invariant state is generated at t = and is annihilated at 
t = T. The four quarks (QQQQ) are spatially fixed in R 3 for 
< t < T. 



assuming that the inter- meson force is subdominant. 



The 4Q Wilson loop is defined by 



W 4Q = -tr (M 1 L 12 M 2 R 12 ). 



(8) 
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Mi, Li.Ri (i=l,2) are given by 



M, 



,Li,Ri = Pexp{ig / 

J A 



dx^A^x)} e SU(3) C .(9) 



Mi.Li.Ri 



Here, L\2 and R\2 are defined by 



jjl a — _ ^abc ^a' b' c' jj)b' j^cc' 



2 • 



£>a' a ^abc^a'b' c' r>bb' r>cc' 



The ground-state 4Q potential V4Q is extracted as 



V4Q = - lim -ln(W 4Q ). 

1 — >-00 1 



(10) 

(11) 

(12) 



In general, the 4Q Wilson loop (W^q) contains excited- 
state contributions, and is expressed as 



(Waq) =Y,C n 



-V n T 



(13) 



where Vo denotes the ground-state 4Q potential V4Q and 
V n (n = 1, 2, 3..) the nth excited-state potential. In prin- 
ciple, V4Q can be obtained from the behavior of (W4Q) at 
the large T region where the ground-state contribution 
becomes dominant. In the practical simulation, however, 
the accurate measurement of V4Q is not easy for large T, 
since (W4Q) decreases exponentially with T. 

To extract the ground-state potential V4q in lattice 
QCD, we adopt the gauge-covariant smearing method 
0- [J EJ to enhance the ground-state compo- 
nent of the 4Q state in the 4Q Wilson loop. The smear- 
ing is known to be a powerful method for the accurate 
measurement of the Q-Q [1 Q and the 3Q potentials 
HL IE 0- an d is expressed as the iterative replacement 
of the spatial link variables Ui(s) (i=l,2,3) by the ob- 
scured link variables Ui(s) S SU(3) C which maximizes 
Retr {UJ(s)Vi(s)} with 



numerical relation Cqq ~ 0.0506a~ 2 obtained from the 
fitting analysis on the on-axis data of the QQ potential 
in lattice QCD at (3 — 6.0 0, @. The gauge configura- 
tions are taken every 500 sweeps after 5000 sweeps using 
the pseudo-heat-bath algorism. We use 300 configura- 
tions for the 4Q potential simulation. For the estimation 
of the statistical error of the lattice data, we adopt the 
jack-knife error estimate. 

On the 16 3 x 32 lattice, we investigate the typical con- 
figuration of 4Q systems as shown in Figs. 6 and 7. In 
Fig. 6, the 4Q system has a planar structure. In Fig. 7, the 
4Q system has a twisted (three-dimension) structure. In 
particular, we analyze in detail the symmetric planar and 
twisted 4Q configurations with d\ — d,2 — = = d, 
although more general asymmetric 4Q configurations 
with various (d±, d2,d$, e^) are also investigated. 

For the 4Q configurations with h < 8, we identify 16 3 
as the spatial size and 32 as the temporal one. On the 
other hand, the calculation for the large-size 4Q config- 
urations with h > 8 is performed by identifying 16 2 x 32 
as the spatial size and 16 as the temporal one. In both 
cases, we use corresponding translational and rotational 
symmetries on lattices for the calculation of (W^q). 
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FIG. 6: A planar configuration of the tetraquark system. 
Q1Q2 is parallel to Q3Q4, and H1H2 is perpendicular to Q1Q2 
and Q3Q4. We call the cases with di = = 1^3 = gLi = d as 
"symmetric planar 4Q configurations". 



V i (s) = aU i (s) + J2T,i U ^ U i( s± ^ U U s + i)} ( 14 ) 

m ± 

with the simplified notation of = Uj(s-j). We here 
adopt the smearing parameter a = 2.3 and the itera- 
tion number N smT = 30, which enhance the ground-state 
component in the 4Q Wilson loop at (3=6.0 in most cases. 
(See the next section.) 



H, 




IV. LATTICE QCD RESULTS FOR THE 4Q 
POTENTIAL 

The lattice QCD simulations are performed with the 
standard plaquette action at (3 = 6.0 on the 16 3 x 32 
lattice at the quenched level. The lattice spacing a is 
estimated as a ~ 0.104fm, which leads to the string ten- 
sion cjqq ~ (427MeV) 2 in the QQ potential, using the 



FIG. 7: A twisted configuration of the tetraquark system. 
Q1Q2 is perpendicular to Q3Q4, and H1H2 is perpendicular 
to Q1Q2 and Q3Q4. We call the cases with di = cfe = ^3 = 
di = d as "symmetric twisted 4Q configurations" . 

For these types of 4Q configurations, we construct the 
4Q Wilson loop W4Q with the junctions locating at Hi 
and H 2 , and calculate the 4Q potential V 4 q from (W40) 
using the smearing method. 
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For the suitable choice of the smearing parameter a 
and the iteration number N smI in Ea. (|14f> . we perform 
some numerical tests with various values of a and iV smr , 
and finally adopt a — 2.3 and N sml - — 30, which are found 
to enhance the ground-state component in the 4Q Wilson 
loop at (3=6.0 in most cases. For the demonstration, 
we show in Fig. 8 a typical example of the ground-state- 
overlap quantity 



C ^(^ 4Q (T)) T +V(W 4Q (T + 1)) 



(15) 



plotted against the iteration number 7V smr at a = 2.3 
for the symmetric planar 4Q configuration with (d, h) = 
(1,5). Here, the ground-state-overlap quantity C indi- 
cates the magnitude of the ground-state component [H,|(| , 
and is found to take a large value close to unity around 
N sjnI = 30. (Note here that the enhancement of the 
ground-state component is the aim of the smearing, and 
hence any approximate optimization is applicable as long 
as the ground-state overlap is enough large.) 




smr 



FIG. 8: A typical example of the ground-state-overlap quan- 
tity Co = <W 4 q(T)> t+1 /(W4q(T+ 1)> T plotted against the 
iteration number N 3ml at a — 2.3 for the symmetric planar 
4Q configuration with (d,h) = (1,5). Co takes a large value 
close to unity around N smr = 30. 

Due to the smearing, the ground-state component 
is largely enhanced in most cases, and therefore the 
4Q Wilson loop (W4Q) composed with the smeared 
link-variable exhibits a single-exponential behavior as 
(W^q) ~ e~ V4QT even for a small value of T. Then, 
for each 4Q configuration, we extract V4Q from the least 
squares fit with the single-exponential form 



(Waq) = Ce 



V iQ T 



(16) 



in the range of T min <T< T max listed in Table I- VI. The 
prefactor C physically means the ground-state overlap, 
and C ~ 1 corresponds to the realization of the quasi- 
ground-state. Here, we choose the fit range of T such 
that the stability of the "effective mass" 



is observed in the range of T m i n <T< T max — 1. 

To see how excited-state contamination is removed in 
this calculation, we show in Fig. 9 several effective-mass 
plots, V(T) v.s. T, for planar and twisted 4Q configura- 
tions at small, intermediate and large distances, respec- 
tively. Owing to the smearing, the effective mass V(T) 
seems to be stable even for small T. To show the qual- 
ity of the single-exponential fit for {W4.Q} as in Ea. dlGII . 
we list the chi-square per degree of freedom, xV-^df, 
for each fit in Table I- VI. For most cases, xV^df takes 
a small value less than unity, and the fitting seems to 
be plausible. (We note that the errors listed in Table 
I- VI are the statistical ones, and there are some sys- 
tematical errors in the lattice QCD calculation. For in- 
stance, when T m ; n = 2 is adopted for large 4Q systems 
as (d, /i)=(3,16), the systematical error originating from 
the fit-range choice seems to be several times larger than 
the statistical error.) 
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FIG. 9: The effective-mass plots, V(T) v.s. T, for several 
4Q configurations at small, intermediate and large distances: 
the symmetric planar and twisted 4Q configurations with 
(d, /i) = ( 1,2), (d, h)=(2,7), and {d, ft)=(3,16). The fit range of 
Tmin < T < T max — 1, which is regarded as the plateau region, 
is indicated by the vertical dashed lines. The solid horizontal 
line denotes the final fit of T4q from the least square fit with 
the single-exponential form. 

In this way, we calculate the tetraquark potential V4Q 
for various 4Q system, i.e., planar, twisted, asymmet- 
ric, and large-size 4Q configurations, respectively. We 
summarize in Table I- VI the lattice QCD results for V4q 
together with the ground-state overlap C, the fit range 
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of T, x 2 /N-df, the minimal flux-tube length L c m ^ for the 
connected 4Q configuration, and the theoretical Ansatze 
V^ 4 h Q and VJqq presented in Sect. II. 

1. Table I and II show V4Q for the symmetric planar 
4Q configurations as shown in Fig. 6 with d\ = d 2 = 
c?3 = di = d. V4Q is shown in terms of d and h. 

2. Table III and IV show V4Q for the symmetric 
twisted 4Q configurations as shown in Fig. 7 with 
d\ = d 2 = d 3 = d 4 = d. V4Q is shown in terms of d 
and h. 

3. Table V shows V4Q for the asymmetric planar 
4Q configurations as shown in Fig. 6 with various 
(^1,(^2)^3,^4) for h=8. 

4. Table VI shows 14q for the asymmetric twisted 
4Q configurations as shown in Fig. 7 with various 
(di, da, da, di) for /i=8. 

Thus, we obtain the tetraquark potential V4Q for about 
200 different patterns of 4Q systems. 

V. DISCUSSIONS 
A. Comparison with theoretical Ansatze 

In this section, we compare the lattice QCD results of 
the 4Q potential V4Q with the theoretical Ansatze pre- 
sented in Sect. II, i.e., the OGE plus multi-Y Ansatz (J5J 
and the two-meson Ansatz (0 . 

For the OGE plus multi-Y Ansatz (JSJ), we set the pa- 
rameters (Aiq, <74q) to be (^3Q, C3Q) in the 3Q potential 
V 3 q in Ref.0, i.e., A 4Q = A 3C , ~ 0.1366, cr 4 Q = <j 3Q ~ 
0.0460a" 2 . Note that there are no adjustable parameters 
except for an irrelevant constant C4Q ~ 1. 2579a -1 . 

For the two- meson Ansatz © , we adopt the lattice re- 
sult for the QQ potential Vqq in Ref.|6j, i.e., ^4qq — 
0.2768, cr Q Q ~ 0.0506a- 2 , C QQ ~ 0.6374a" 1 . Then, 
there are no adjustable parameters also for the two-meson 
Ansatz. 

We demonstrate the two Ansatze for the symmetric 
planar 4Q configurations as shown in Fig. 6 with d\ = 
d 2 = d-3 = d 4 = d. In this case, the OGE Coulomb plus 
multi-Y Ansatz for the connected 4Q system reads 

+ cr4Q_L m i n + C4Q. (18) 

In the case of h > -^d, the lowest connected 4Q system 
takes the double- Y-shaped flux-tube, and the minimal 
value of the total flux-tube length is expressed as 

L min = h+ 2y/M. (19) 

In the case of h < S^d, the lowest connected 4Q system 
takes the X-shaped flux-tube with L m i n = 2\/h 2 + 4d 2 , 



although it must be unstable against the decay into two 
mesons. On the other hand, the two-meson Ansatz reads 

V 2Q Q(h) = 2xV QQ (h), (20) 

which is independent of d. 

We show in Fig. 10(a) the lattice QCD results of the 
4Q pot ential V4Q for symmetric planar 4Q configurations 
Hi IM E3 with d = 1 - 4. The symbols denote lattice 
QCD results. The curves describe the theoretical form: 
the solid curve denotes the OGE plus multi-Y Ansatz JSJ), 
and the dashed-dotted curve the two- meson Ansatz ijfJJl. 

For large value of h compared with d, the lattice data 
seem to coincide with the OGE Coulomb plus multi-Y 
Ansatz [12, ^J. On the other hand, for small h, the lat- 
tice data tend to agree with the two-meson Ansatz and 
seem independent of d j38l . These tendencies were 
also observed in a recent lattice work by other group 
|42|. This would correspond to the transition from the 
connected 4Q state into the two-meson state as h de- 
creases, as will be discussed in the next subsection. 

Here, we comment on the transition in terms of the 
ground-state overlap C listed in Table I. For large h, the 
ground-state overlap C is almost unity, which implies 
the realization of the quasi-ground-state in the present 
calculation with the smeared 4Q Wilson loop based on 
the connected 4Q configuration. For small h, however, 
C tends to be small, and hence, for accurate measure- 
ments, we have to take relatively large values of T as 
the fit range. This would indicate that the ground-state 
configuration is largely different from the connected 4Q 
configuration for small h. (In other words, it may be 
nontrivial to obtain the result indicating the two-meson 
state for small h from the 4Q Wilson loop based on the 
connected 4Q configuration.) 

Next, we investigate the twisted 4Q configuration [3^, 
l4lj as shown in Fig. 7. We show in Fig. 10(b) the lattice 
QCD results of the 4Q potential for symmetric twisted 
4Q configurations with d = 1 — 4. The symbols denote 
lattice QCD results for each d, and the curves describe 
the theoretical form of the OGE plus multi-Y Ansatz. 

The lattice data seem to agree with the OGE plus 
multi-Y Ansatz in the wide region of h [3^, ^J. In 
the twisted 4Q configuration, the distance between the 
nearest quark and antiquark cannot take a smaller value 
than the "inter-diquark distance" h, and therefore T4 4 q 
is smaller than V 2 qq in most cases except for extreme 
configurations as d > h. (See Fig. 10(b).) Then, different 
from the planar case, it is not easy to make the tran- 
sition from the connected 4Q state into the two-meson 
state for the twisted case. Also from the lattice data, the 
ground-state overlap C is found to be almost unity for 
all twisted 4Q configurations, which indicates that the 
ground state resembles a connected 4Q state. In other 
words, the twisted 4Q configuration seems to be rather 
stable against the transition into the two QQ mesons, 
which may indicate a stability of the "twisted structure" 
or the "tetrahedral structure" of the 4Q system. 
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FIG. 10: The tetraquark (4Q) potential T4q: (a) for sym- 
metric planar 4Q configurations as shown in Fig.6; (b) for 
symmetric twisted 4Q configurations as shown in Fig. 7. The 
symbols denote the lattice QCD results. The curves describe 
the theoretical form: the solid curve denotes the OGE plus 
multi-Y Ansatz, and the dashed-dotted curve the two-meson 
Ansatz. 



We also investigate more general asymmetric 4Q con- 
figurations with various (di, da> <^3» ^4) for both planar 
and twisted cases, as shown in Table V and VI. Also 
for the asymmetric planar and twisted 4Q configurations, 
V4Q seems to be well described with the OGE plus multi- 
Y Ansatz in the case of h > -^di (i = 1,2,3,4). Note 
here that some 4Q configurations are physically equiva- 
lent, e.g., the planar cases with (dx, di 1 ^3, di)=(l,l,l,2) 
and (0,2,1,2), although the corresponding smeared 4Q 
Wilson loops are different. For such cases, the lattice 
QCD results are found to be almost the same. In fact, 
the extracted lattice results are almost independent of 
the way how the 4Q Wilson loop is constructed, as long 
as the spatial locations of the static four quarks are the 
same. This indicates that the ground-state contribution 
is properly extracted in the present calculation. 

As the conclusion, the OGE plus multi-Y Ansatz well 
describes the 4Q potential V4Q, when QQ and QQ are 
well separated, e.g., the "inter-diquark distance" h is 
large in comparison with the "diquark size" d. On the 
other hand, when the nearest quark and antiquark pair is 
spatially close, the system is described as a "two-meson" 
state. 



Together with the previous studies |alalM|3J,l38, 39, 
E3>E| for the inter-quark potentials in lattice QCD, we 
have found the universality of the string tension as 

(Tqq ~ 03Q ~ CT 4Q ~ CJ5Q - (420MeV) 2 , (21) 

and the OGE result for the Coulomb coefficient as 



^QQ - 2A 3Q 



2A 



4Q 



2A 



5Q 



0.27. 



(22) 



In particular, these lattice QCD studies [36|, |3JJ |38j, |3tj, 

l40l Hlj indicate a fairly good agreement among CT3Q , <t 4 q 
and 05Q , which seem to be slightly smaller than ctqq . (As 
an interesting possibility, the numerical similarity among 
03Q, CT4Q and 05Q may reflect the similar structure of 
the Y-type flux-tube in the multi-quark systems.) The 
universality of the string tension observed in our lattice 

QCD studies HHIHilSlilliiliS SeemS t0 be 

consistent with the hyp othetical flux-tube picture |49Ll50l 
IFll I52I l53l l54l l55l l56j or the strong-coupling expansion 
scheme [5(J, HH , although strong-coupling QCD does not 
have a continuum limit and is far from real QCD. As 
for the irrelevant constant, Cqq, C3Q, C4Q and C5Q are 
non-scaling unphysical quantities appearing in the lattice 
regularization, and we find an approximate relation as 



C 



3Q 



a 



4Q 



a 



5Q 



in lattice QCD 



0.32a" 



(23) 



B. The quark confinement force in 4Q systems 

While the short-distance OGE Coulomb force can be 
understood with perturbative QCD, the long-distance 
confinement force is a typical nonperturbativc quantity 
and highly nontrivial particularly for multi-quark sys- 
tems. To specify the long-distance property of V4Q is 
important to clarify the confinement mechanism from a 
wide viewpoint including multi-quarks, and it also leads 
to a proper quark-model Hamiltonian to describe multi- 
quark systems. Therefore, we perform a further analysis 
for the long-distance force in 4Q systems. 

To clarify the long-distance force in the 4Q system, we 
plot V4Q against L^-J for planar and twisted 4Q config- 
urations in Figs. 11 (a) and (b), respectively. Here, 1/^2 
denotes the minimal flux-tube length for the connected 
4Q system. In both planar and twisted cases, for large 
-^mS' ^4Q approaches a linearly arising function of 

To single out the long-distance confinement force, we 
consider the 4Q potential subtracted by the Coulomb 
part. Here, we subtract the OGE Coulomb part V^q 
of V C 4Q in Eq.JSJ) for the connected 4Q system, with the 
Coulomb coefficient A4Q fixed to be A$q in the 3Q poten- 
tial y 3 Q in Rcf.0. We plot V iQ - V^ 1 against L c *9 for 
planar and twisted 4Q configurations in Figs. 12 (a) and 

^Cou" 
c4Q 

C4Q except for a small h region, 



(b) , respectively. For the planar 4Q system, V4Q — K? oul 
approaches ct 4 qL^ 
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FIG. 11: The 4Q potential V 4 q plotted against L™J: (a) for 
planar 4Q configurations; (b) for twisted 4Q configurations, 
-^min denotes the minimal flux-tube length for the connected 
4Q configuration. We plot all the lattice QCD data of V4Q 
for the symmetric case with d = 1 (solid circles), d = 2 (solid 
squares), d — 3 (open triangles) and d = 4 (crosses) together 
with the asymmetric case (open diamonds). 



FIG. 12: The 4Q potential subtracted by the OGE Coulomb 
part of the connected 4Q system, T4q — V C 4q u1 , plotted against 
_Lmin : (a) for planar 4Q configurations; (b) for twisted 4Q 
configurations. We plot all the lattice QCD data including 
asymmetric cases. The meaning of the symbols is the same 
as in Fig.ll. The solid line denotes o-^qL^® + C4Q. 



where the flip-flop into a two-meson state can take place. 
For the twisted 4Q system, we observe remarkable agree- 
ment between the lattice QCD data of V4Q — V^q" 1 and 

C4qL£J2 + ^4Q f° r the wide region of L m ; n , which cor- 
responds that the flip-flop into the two-meson state does 
not occur in most twisted 4Q configurations. 

Thus, the confinement potential in the 4Q system as 
shown in Fig. 2 is proportional to L mul , which indicates 
that the quark confinement force is genuinely 4-body and 
the flux tube is multi-Y shaped. 



C. The flip-flop, the level crossing, and absence of 
the color van der Waals force 

Finally, we investigate the flip-flop between the con- 
nected 4Q state and the two-meson state. Since the flux- 
tube changes its shape so as to have the minimal length, 
the multi-Y type flux-tube is expected to change into a 
two-meson state for small h. 

This type of the flip-flop is physically important for the 
properties of 4Q states especially for their decay process 
into two mesons. Note also that, in the flux-tube picture, 
the meson-meson reaction is described by the flux-tube 



recombination between the two mesons, and this process 
can be realized through the two successive flip-flops be- 
tween the two-meson state and the connected 4Q state. 
Therefore, this type of the flip-flop is important also for 
the reaction mechanism between two mesons. 

As a clear signal of the flip-flop, we again show the 4Q 
potential V4Q for the symmetric planar 4Q configuration 
with d=l,2,3 separately in Fig. 13. The solid curve de- 
notes V C 4Q for the OGE plus multi-Y Ansatz, and the 
dashed-dotted curve V^qq = 2Vqq{K) for the two- meson 
Ansatz. For large h, V4Q coincides with the energy V C 4Q 
of the connected 4Q system. For small h, V4Q coincides 
with the energy 2Vqq of the "two-meson" system com- 
posed of two flux-tubes. In the intermediate region of 
h, one can observe the cross over from one Ansatz to 
another. 

Thus, in these particular cases, we can observe a clear 
evidence of the flip-flop as 

V4 Q (d,h)~mm(V C 4 Q (d,h),2V QQ (h)), (24) 

which indicates the transition between the connected 
4Q state and the two-meson state around the level- 
crossing point where these two systems are degenerate 
as V C 4Q(d,h) = 2Vqq(/i). This result also supports the 
flux-tube picture even for the 4Q system. 
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The present lattice QCD results on the "flip-flop" leads 
to infrared screening and disappearance of the long-range 
color interactions, i.e., the confining force and the OGE 
Coulomb force, between (anti)quarks belonging to dif- 
ferent "mesons". This physically results in the absence 
of the tree-level color van der Waals force between two 
mesons El El El. 
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FIG. 13: The typical lattice QCD results for the flip-flop be- 
tween the connected 4Q state and the two-meson state for the 
symmetric planar 4Q configuration with fixed d. The symbols 
denote lattice QCD results. The curves describe the theo- 
retical form: the solid curve denotes the OGE plus multi-Y 
Ansatz, and the dashed-dotted curve the two-meson Ansatz. 



at the quenched level. For about 200 different patterns 
of 4Q systems, we have extracted V4Q from the 4Q Wil- 
son loop in 300 gauge configurations, with the smear- 
ing method to enhance the ground-state component. We 
have calculated V4Q for planar, twisted, asymmetric, and 
large-size 4Q configurations, respectively. The calcula- 
tion for large-size 4Q configurations has been done by 
identifying 16 2 x 32 as the spatial size and 16 as the tem- 
poral one, and the long-distance confinement force has 
been particularly analyzed in terms of the flux-tube pic- 
ture. 

When QQ and QQ are well separated, V4Q is found 
to be expressed as the sum of the one-gluon-exchange 
Coulomb term and multi-Y type linear term based on 
the flux-tube picture. In this case, all the four quarks 
are linked by the connected double- Y-shaped flux-tube, 
where the total flux-tube length is minimized. On the 
other hand, when the nearest quark and antiquark pair is 
spatially close, the system is described as a "two-meson" 
state rather than the connected 4Q state. 

We have observed a flux-tube recombination called as 
"flip-flop" between the connected 4Q state and the "two- 
meson" state around the level-crossing point. This "flip- 
flop" leads to infrared screening of the long-range color 
interactions between (anti)quarks belonging to different 
mesons, and results in the absence of the tree-level color 
van der Waals force between two mesons. 

As a next step, it is interesting to investigate the tran- 
sition in terms of the level crossing between the connected 
4Q state and the two-meson state through the diagonal- 
ization of the correlation matrix with various 4Q states 
[34| . Through the investigation of the excited-state lev- 
els of the 4Q system, a realistic picture for the reaction 
mechanism between two mesons may be obtained. 

The dynamical quark effect for the flux-tube picture 
and the flip-flop would be also an interesting subject. In 
this context, the string-breaking effect may cause more 
complicated variation of the transition between single 
QQ meson and the multi-quark system. 

In any case, recent lattice QCD studies begin to shed 
light on the realistic picture in hadron physics and to 
reveal even the properties of the multi-quark system. 
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APPENDIX A: OGE COULOMB TERMS IN THE 
4Q POTENTIAL 

In the Appendix, we briefly show the derivation of the 
OGE Coulomb terms in the 4Q potential in Eq.JSJ) by 
calculating (T|TfT?|T) for the tetraqurak (4Q) state |T). 

In the quark picture, the 4Q state |T) corresponding 
to Fig. 2 is expressed as 



IT) = |(QiQ 3 )s(Q3Q4)3)i 

r abc r ab' c 



2V3 



IQ1Q2Q3 Q4 



(Al) 



where the indices denote the SU(3) color indices of 
(anti)quarks. Here, |T) is normalized as 



(T|T) 



abc ab c def de f 
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1. 



(A2) 



The color matrix factor T^T" in the OGE process can 
be expressed with the Casimir operator C2(R) as 



After some calculation, one finds 



|ii) = |(QiQ 3 )i(Q2Q 4 )i) - ^QIQlQlQ 



^&>, (A6) 



|88) = |(QiQ 3 ) 8 (Q2Q4: 



8 1 



= ^{|Q;oSqSQ2> - 5lQ?QlQIQ|», (A7) 

which satisfy the orthonormal condition, 

(11|11) = (88|88) = 1, <11|88) = 0. (A8) 

Using Egs- lHH , C^l 1 and 113; c i and c % can be ob- 
tained as 



Ci = <11|T) 



1 
1 

6V3 



,abc _ab' c 



(Q1QIQ3Q4IQ1Q2Q3Q4 



^abc^abc 



1 

V5' 



(A9) 



rparpa 

i 3 



1 

2 
1 

2 

-C 2 (Ri+j) - -, 



(Tt) 2 

{C2(Ri+j) — C 2 (.Rj) - 



(A3) 



where Ri+j denotes the total SU(3) color representation 
of the (i + j) system. Here, C 2 (i?i) — C2(Rj) — 3 has 
been used for each (anti)quark belonging to 3 (3). 

In this 4Q system, the two quarks, Qi and Q2, form the 
3 representation, i.e., C 2 (i?i +2 ) = C 2 (3) = |, and then 
one gets (T|T"T^|T) = -§. This type of the Coulomb 
coefficient between two quarks is the same as that in the 
3Q system. Similarly, one gets (T|T 3 a T 4 a |T) = -§ for the 
two antiquarks, Q3 and Q4. 

Next, we consider the Coulomb interaction between 
the quark and the antiquark. Owing to the symmetry, 
we only have to investigate the interaction between Qi 
and Q3. To this end, we rewrite the 4Q state |T) in terms 
of the irreducible representation fo the Q1+Q3 system. 
Since Qi and Q 3 can form the singlet (1) or the octet (8) 
representation, the 4Q state |T) can be rewritten as 



|T) = Ci|(QiQ 3 )i(Q 2 Q4)i> 

+ C 8 |(QlQ3)8(Q 2 Q4)8>l 

with appropriate constants C\ and Cg satisfying 



(A4) 



C 8 = (88|T) 



- (Q1Q2Q3Q4IQ1Q2Q3 Q4 



^abc^ab c / 



1 



12-s/6 
1 



(QiQ 2 Q3Q4lQiQ 2 Q3 Q4 



^abc^acb 



abc abc 



= -\/g- (MO) 



4V6 ' nV6 
Then, using Eq.(£j), C 2 (l) = and C 2 (8) = 3, we get 

(T|TfT 3 °|T) = IC^i^lllTfl^Ill) + |(7 8 | 2 (88|rfr^|88> 
= ICxI^lc^l) - |} + |C 8 | 2 {^C 2 (8) - |} 



(AH) 



In this way, for the Coulomb interaction in the 4Q 
system as shown in Fig. 3, we obtain 



(T\T?T%\T) = (T|T 3 a T 4 a |T) = --, (A12) 
(T\TiTs\T) = (T|T"T4 |T) = (T|T 2 a T 3 Q |T) 

= (T|T 2 a T 4 a |T) = -i (A13) 



Idl 



|C 8 | 2 = 1. 



(A5) and derive the Coulomb terms in Eq.0. 



[1] M. Creutz, Phys. Rev. Lett. 43, 553 (1979), Erratum- edition (World Scientific, 1997) and references therein. 

ibid. 43, 890 (1979); Phys. Rev. D21, 2308 (1980). [3] APE Collaboration (M. Albanese et al), Phys. Lett. 

[2] For instance, H. J. Rothe, "Lattice Gauge Theories", 2nd B192, 163 (1987). 



12 



[4] G.S. Bali and K. Schilling, Phys. Rev. D47, 661 (1993). 

[5] T. T. Takahashi, H. Matsufuru, Y. Nemoto, and H. Sug- 
anuma, Phys. Rev. Lett. 86, 18 (2001). 

[6] T. T. Takahashi, H. Suganuma, Y. Nemoto, and H. Mat- 
sufuru, Phys. Rev. D 65, 114509 (2002). 

[7] T. T. Takahashi and H. Suganuma, Phys. Rev. Lett. 90, 
182001 (2003); Phys. Rev. D70, 074506 (2004). 

[8] H. Suganuma, H. Ichie, and T.T. Takahashi, Proc. of Int. 
Conf. on "Color Confinement and Hadrons in Quantum 
Chromodynamics" , July 2003, RIKEN, Japan, edited by 
H. Suganuma et al, (World Scientific, 2004) 249-261. 

[9] LEPS Collaboration (T. Nakano et al), Phys. Rev. Lett. 
91, 012002 (2003). 
[10] DIANA Collaboration (V. V. Narmin et al), Phys. Atom. 

Nucl. 66, 1715 (2003). 
[11] CLAS Collaboration (S. Stephanian et al.), Phys. Rev. 

Lett. 91, 252001 (2003). 
[12] SAPHIR Collaboration (J. Barth et al), Phys. Lett. 

B572, 127 (2003). 
[13] HI Collaboration (A. Aktas et al), Phys. Lett. B588, 17 
(2004). 

[14] NA49 Collaboration (C. Alt et al), Phys. Rev. Lett. 92, 
042003 (2004). 

[15] D. Diakonov, V. Petrov and M. Polyakov, Z. Phys. A359, 
305 (1997). 

[16] For a recent review, S. L. Zhu, Int. J. Mod. Phys. A19, 
3439 (2004) and references therein. 

[17] BES Collaboration (J.Z. Bai et al), Phys. Rev. D70, 
012004 (2004); HERA-B Collaboration (I. Abt et al), 
Phys. Rev. Lett 93, 212003 (2004); SPHINX Collabo- 
ration (Y.-M. Antipov et al), Eur. Phys. J. A21, 455 
(2004); HyperCP Collaboration (M.J. Longo et al), 
Phys. Rev. D70, 111101 (2004). 

[18] Belle Collaboration (S. K. Choi et al), Phys. Rev. Lett. 
91, 262001 (2003). 

[19] CDF II Collaboration (D. Acosta et al), Phys. Rev. Lett. 
93, 072001 (2004). 

[20] DO Collaboration (V. M. Abazov et al), Phys. Rev. Lett. 
93, 162002 (2004). 

[21] BABAR Collaboration (B. Aubert et al), Phys. Rev. 
Lett. 93, 041801 (2004). 

[22] BABAR Collaboration (B. Aubert et al), Phys. Rev. 
Lett. 90, 242001 (2003). 

[23] Belle Collaboration (P. Krokovny et al), Phys. Rev. Lett. 
91, 262002 (2003). 

[24] F. E. Close and S. Godfrey, Phys. Lett. B574, 210 (2003). 

[25] F. E. Close, P. R. Page, Phys. Lett. B578, 119 (2004). 

[26] N. A. Tornqvist, "Comment on the Narrow Charmo- 
nium State of Belle at 3871.8-MeV as a Deuson", 
hep-ph/0308277 (2003). 

[27] S. Pakvasa and M. Suzuki, Phys. Lett. B579, 67 (2004). 

[28] C. Y. Wong, Phys. Rev. C69, 055202 (2004). 

[29] E. Braaten and M. Kusunoki, Phys. Rev. D69, 074005 
(2004). 

[30] T. Barnes and S. Godfrey, Phys. Rev. D69, 054008 
(2004). 

[31] E. J. Eichten, K. Lane and C. Quigg, Phys. Rev. D69, 
094019 (2004). 

[32] E. S. Swanson, Phys. Lett. B588, 189 (2004); ibid. B598, 
197 (2004). 

[33] R.L. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 
(2003). 

[34] A. M. Green, J. Lukkarinen, P. Pennanen, and 
C. Michael, Phys. Rev. D 53, 261 (1996). 



[35] P. Pennanen, A. M. Green, and C. Michael, Phys. Rev. 

D 59, 014504 (1999). 
[36] F. Okiharu, H. Suganuma, and T. T. Takahashi, "First 

study for the pentaquark potential in SU(3) lattice 

QCD", hep-lat/0407001, Phys. Rev. Lett. 94, 192001 

(2005). 

[37] H. Suganuma, T. T. Takahashi, F. Okiharu, and H. Ichie, 
Proc. of "QCD Down Under", March 2004, Adelaide, 
Nucl. Phys. B (Proc. Suppl.) 141, 92 (2005). 

[38] H. Suganuma, T. T. Takahashi, F. Okiharu, and H. Ichie, 
Proc. of "Quark Confinement and Hadron Spectrum VI", 
Sept. 2004, Italy, edited by N. Brambilla et al, AIP Conf. 
Proc. CP756, 123 (2005). 

[39] H. Suganuma, F. Okiharu, T. T. Takahashi, and H. Ichie, 
Nucl. Phys. A755, 399 (2005). 

[40] H. Suganuma, H. Ichie, F. Okiharu, and T. T. Takahashi, 
Proc. of Int. Workshop on "Pentaquark04" , July 2004, 
SPring-8, Japan (World Scientific, 2005) 414. 

[41] F. Okiharu, H. Suganuma, and T. T. Takahashi, Proc. of 
Int. Workshop on "Pentaquark04" , July 2004, SPring-8, 
Japan (World Scientific, 2005) 339. 

[42] C. Alexandrou, and G. Koutsou, Phys. Rev. D71, 014504 
(2005). 

[43] J. Weinstein and N. Isgur, Phys. Rev. Lett. 48, 659 
(1982). 

[44] Fl. Stancu and D. O. Riska, Phys. Lett. B575, 242 
(2003). 

[45] Y. Kanada-Enyo, O. Morimatsu, and T. Nishikawa, 

Phys. Rev. C71, 045202 (2005); Proc. of Int. Workshop 

on "Pentaquark04" , July 2004, SPring-8, Japan, (World 

Scientific, 2005) 239. 
[46] R. L. Jaffe, Phys. Rev D 15, 267 (1977). 
[47] D. Black, A. H. Fariborz, F. Sannino, and J. Schechter, 

Phys. Rev. D 59, 074026 (1999). 
[48] M. Alford and R. L. Jaffe, Nucl. Phys. B578, 367 (2000). 
[49] Y. Nambu, Phys. Rev. D10, 4262 (1974). 
[50] J. Kogut and L. Susskind, Phys. Rev. Dll, 395 (1975). 
[51] J. Carlson, J. Kogut, and V. Pandharipande, Phys. Rev. 

D27, 233 (1983); D28, 2807 (1983). 
[52] A. Casher, H. Neuberger, and S. Nussinov, Phys. Rev. 

D20, 179 (1979). 
[53] N. Isgur and J. Paton, Phys. Lett. B124, 247 (1983); 

Phys. Rev. DB31, 2910 (1985). 
[54] A. Patel, Nucl. Phys. B243, 411 (1984); Phys. Lett. 

B139, 394 (1984). 
[55] M. Oka, Phys. Rev. D 31, 2274 (1985); M. Oka and 

C.J. Horowitz, Phys. Rev. D31, 2773 (1985). 
[56] F. Lenz, J.T. Londergan, E.J. Moniz, R. Rosenfelder, 

M. Stingl, and K. Yazaki, Annals Phys. 170, 65 (1986). 
[57] H. Ichie, V. Bornyakov, T. Streuer, and G. Schierholz, 

Nucl. Phys. A721, 899 (2003); Nucl. Phys. B (Proc. 

Suppl.) 119, 751 (2003). 
[58] V.G. Bornyakov, H. Ichie, Y. Mori, D. Pleiter, M.I. Po- 

likarpov, G. Schierholz, T. Streuer, H. Stiiben, and 

T. Suzuki, Phys. Rev. D70, 054506 (2004). 
[59] P.O. Bowman and A. P. Szczepaniak, Phys. Rev. D70, 

016002 (2004). 

[60] D.S. Kuzmenko and Yu. A. Simonov, Phys. Atom. Nucl. 

66, 950 (2003). 
[61] J.M. Cornwall, Phys. Rev. D69, 065013 (2004). 
[62] T. Appelquist and W. Fischler, Phys. Lett. B77, 405 

(1978). 

[63] M.B. Gavela, A. Le Yaouanc, L. Oliver, O. Pene, 
J.C. Raynal, and S. Sood, Phys. Lett. B82, 431 (1979). 



13 

[64] G. Feinberg and J. Sucher, Phys. Rev. D20, 1717 (1979). 



14 



TABLE I: A part of lattice QCD results of the 4Q potential V4Q for the planar 4Q configuration as shown in Fig. 6 with 
di — d,2 — do, = d,4 = d. The symmetric planar 4Q systems are labeled by (d, h). We list also the ground-state overlap C, the fit 
range of T, x 2 /Nuf, the minimal flux-tube length for the connected 4Q configuration, and the theoretical Ansatze. V^q 
denotes the OGE Coulomb plus multi-Y Ansatz J5J with (A4q,(T4q) fixed to be (A^Q,a^Q) in V3Q in Ref.|jj. V^qq denotes the 
two-meson Ansatz JSJ with Vqq in Ref.|6J. All the quantities are measured in the lattice unit at f3 — 6.0, i.e., a ~ 0.104fm. 



(d,h) 


V4Q 


C 


i- ruin '-/ niax 


x7-Wdf 


£C4Q 
min 


V th 

c4Q 


V th - 

2QQ 


(1,1) 

V- 1 -' ^ } 


0.8590(45) 


0.4970(135) 


6-11 


0.140 


4.47 


1.1293 


0.8224 


(1,2) 


1.2124(31) 


0.9049(111) 


4-9 


0.710 


5.46 


1.2560 


1.2004 


(1,3) 


1.3218(17) 


0.9572(48) 


3-9 


0.195 


6.46 


1.3352 


1.3939 


(1,4) 


1.3938(11) 


0.9627(19) 


2-8 


0.755 


7.46 


1.3999 


1.5412 


(1 5) 


1.4531(13) 


0.9556(22) 


2-8 


0.633 


8.46 


1.4579 


1.6701 


(1,6) 


1.5051(33) 


0.9354(95) 


3-9 


0.072 


9.46 


1.5123 


1.7897 


(1 7) 
V- 1 -, 1 1 


1.5644(19) 


0.9434(32) 


2-7 


0.336 


10.46 


1.5643 


1.9041 


(1,8) 


1.6184(20) 


0.9375(36) 


2-8 


0.124 


11.46 


1.6150 


2.0152 


(1,9) 


1.6706(22) 


0.9295(35) 


2-4 


0.119 


12.46 


1.6646 


2.1241 


(1,10) 


1.7269(28) 


0.9286(48) 


2-7 


0.184 


13.46 


1.7136 


2.2314 


(1,11) 


1.7774(30) 


0.9177(48) 


2-4 


0.001 


14.46 


1.7620 


2.3377 


(1,12) 


1.8234(37) 


0.8983(62) 


2-7 


1.094 


15.46 


1.8100 


2.4431 


(1 13) 


1.8797(44) 


0.8979(74) 


2-7 


0.351 


16.46 


1.8577 


2.5478 


(1 14) 


1.9285(49) 


0.8838(78) 


2-6 


0.899 


17.46 


1.9052 


2.6521 


(1,15) 


1.9787(55) 


0.8722(91) 


2-7 


0.361 


18.46 


1.9525 


2.7559 


(1,16) 


2.0365(62) 


0.8741(101) 


2-6 


0.702 


19.46 


1.9996 


2.8594 


(2,1) 


0.8263(46) 


0.3324(93) 


6-13 


0.624 


8.25 


1.3992 


0.8224 


(2,2) 


1.2633(183) 


0.5646(514) 


5-9 


0.727 


8.94 


1.5022 


1.2004 


(2,3) 


1.4592(111) 


0.7457(328) 


4-7 


0.119 


9.93 


1.5734 


1.3939 


(2,4) 


1.5950(17) 


0.9223(29) 


2-7 


0.643 


10.93 


1.6340 


1.5412 


(2,5) 


1.6619(21) 


0.9286(36) 


2-6 


0.779 


11.93 


1.6896 


1.6701 


(2,6) 


1.7215(23) 


0.9285(41) 


2-7 


0.653 


12.93 


1.7426 


1.7897 


(2,7) 


1.7791(29) 


0.9283(51) 


2-6 


0.089 


13.93 


1.7938 


1.9041 


(2,8) 


1.8279(31) 


0.9121(54) 


2-6 


0.631 


14.93 


1.8439 


2.0152 


(2,9) 


1.8827(34) 


0.9101(57) 


2-7 


0.111 


15.93 


1.8932 


2.1241 


(2,10) 


1.9366(38) 


0.9049(64) 


2-6 


0.108 


16.93 


1.9419 


2.2314 


(2,11) 


1.9917(45) 


0.9026(73) 


2-5 


0.785 


17.93 


1.9902 


2.3377 


(2,12) 


2.0419(49) 


0.8911(82) 


2-5 


0.215 


18.93 


2.0381 


2.4431 


(2,13) 


2.0954(55) 


0.8856(90) 


2-6 


0.445 


19.93 


2.0857 


2.5478 


(2,14) 


2.1387(61) 


0.8619(101) 


2-5 


1.272 


20.93 


2.1331 


2.6521 


(2,15) 


2.1990(67) 


0.8683(111) 


2-5 


1.103 


21.93 


2.1804 


2.7559 


(2,16) 


2.2482(78) 


0.8567(122) 


2-5 


0.687 


22.93 


2.2274 


2.8594 
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TABLE II: A part of lattice QCD results of the 4Q potential Vaq for the symmetric planar 4Q configuration labeled by (d, h). 
The notations are the same as in Table I. 



(d, h) 


ViQ 


C 




X 2 /Ndf 


L c4Q 
mm 


V th 


V th - 

2QQ 


(3 11 


8281f20l 


3242(331 


5-12 


0.502 


12.17 


1.6129 


0.8224 




1.2143(228) 


0.3390(387) 


5-9 


0.237 


12.65 


1.7043 


1.2004 


(3,3) 


1.5031(168) 


0.5540(370) 


4-7 


2.112 


13.42 


1.7636 


1.3939 


(3 41 


1 6992(681 


9 7751 (1581 


3-7 


0.654 


14.39 


1.8213 


1.5412 


(3,5) 


1 81 85(31 1 


f) 8797(531 


2-6 


0.446 


15.39 


1.8756 


1.6701 


(3,6) 


1.8896(35) 


0.8922(60) 


2-6 


0.299 


16.39 


1.9275 


1.7897 


("3 7) 


1.9459(37) 


0.8849(63) 


2-4 


1.708 


17.39 


1.9780 


1.9041 


(3,8) 


2 nm 6(451 


f) 8898(761 


2-5 


0.933 


18.39 


2.0276 


2.0152 


(3,9) 


2.0544(45) 


0.8739(73) 


2-4 


0.003 


19.39 


2.0766 


2.1241 


(3 10) 


2 1123(521 


D 8748(851 


2-6 


0.116 


20.39 


2.1250 


2.2314 


(3,11) 


2 1658(571 


8698(921 


2-7 


0.241 


21.39 


2.1731 


2.3377 


(3,12) 


2 2234(651 


9 8720(1051 


2-5 


0.373 


22.39 


2 2208 


2.4431 


(3,13) 


2 2698(731 


8548(1 181 


2-5 


0.168 


23.39 


2 2683 


2.5478 


(3,14) 


2 3192(881 


8418(1431 


2-5 


0.417 


24.39 


2.3157 


2.6521 


(3,15) 


2 3843(941 


8564(1541 


2-6 


0.653 


25.39 


2.3628 


2.7559 


(3,16) 


2.4393(112) 


8530(1851 


2-5 


0.469 


26.39 


2.4098 


2.8594 


(4 11 


0.8228(21) 


0.3069(34) 


5-12 


0.586 


16.12 


1.8119 


0.8224 


(4 21 


1 2510(771 


3541 (1 1 01 


4-9 


0.528 


16.49 


1.8975 


1.2004 


(4 31 


1 4643(2471 


3390(3361 


4-7 


0.247 


17.09 


1.9482 


1.3939 


(4 4) 


1 7452(1 121 


5783(1941 


3-7 


0.910 


17.89 


1.9972 


1.5412 


(4 51 


1.9147(159) 


0.7137(337) 


3-5 


1.694 


18.86 


2.0493 


1.6701 


(4 61 


2 0168(2301 


7591(5181 


3-7 


0.425 


19.86 


2.1007 


1.7897 


(4 71 


2 1167(491 


8435(841 


2-4 


0.041 


20.86 


2.1507 


1.9041 


(4 81 


2.1712(56) 


0.8341(89) 


2-4 


0.436 


21.86 


2.2000 


2.0152 


(4,9) 


2.2444(62) 


0.8611(102) 


2-6 


0.363 


22.86 


2.2486 


2.1241 


(4,10) 


2.2970(84) 


0.8505(138) 


2-5 


1.626 


23.86 


2.2968 


2.2314 


(4,H) 


2.3430(84) 


0.8332(133) 


2-6 


0.514 


24.86 


2.3447 


2.3377 


(4,12) 


2.3976(96) 


0.8309(149) 


2-5 


0.544 


25.86 


2.3923 


2.4431 


(4,13) 


2.4503(109) 


0.8225(168) 


2-4 


0.030 


26.86 


2.4397 


2.5478 


(4,14) 


2.5043(123) 


0.8178(194) 


2-4 


0.070 


27.86 


2.4869 


2.6521 


(4,16) 


2.6037(177) 


0.7952(276) 


2-4 


2.015 


29.86 


2.5809 


2.8594 
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TABLE III: A part of lattice QCD results of the 4Q potential Vaq for the twisted 4Q configuration as shown in Fig. 7 with 
di = d 2 — d 3 = d 4 = d. The symmetric twisted 4Q systems are labeled by (d, h). The notations are the same as in Table I. 



(d, h) 


Vac^ 
v 4Q 


c 


T ■ — T 

J- mm max 


v 2 /AT 
X / iV DF 


L c4Q 
min 


v c4Q 


V th - 

2QQ 




i.i ( /y(Uo) 


n nsnc ( i i \ 

u.yoyo^ii J 


Z-5 


U.oyo 


4.4/ 


1 1 coo 

l.loyo 


1 1 one 
l.loUo 




l.ZO / / ^UO j 


u.yoo 


z-0 


n nno 
u.uuy 


0.40 


1 9fi1 1 
l.zOll 


1 9Qfi7 

i.zyo / 




l.ooll^Uo ) 


u.yo / 0^10 ) 


z-o 


n /i fi^ 

U.400 


fi /ifi 
0.40 


1 QQfi9 
l.OOOZ 


1 AA'iK 
1.4400 


(M) 


i QGfifi/'l 1 \ 
l.oyoU^ll ) 


u.yo4z(^zu j 


z-O 


n fi7fi 
U.O / 


7 /Ifi 
/ .40 


1 a nno 
1.4UUZ 


1 K7Q7 

1.0 / / 




1 /] EQO/1 Q\ 

1.40oz^lo J 


u.yo4o^zo j 


z-o 


n r^/i fi 
U.040 


O.40 


1 A Kfifl 

1.40oU 


1 fiO/1 1 

i.oy4i 




i.oiuu^oo ) 


u.y4y / i^yy j 


O-O 


n 991 

U.zZl 


Q /Ifi 

y.40 


1 K1 9Q 

i.oizo 


l.oUoo 


(M) 


1 C^l /I Q\ 

l.OOOl^lo ) 


f) O/l 70^ Q0^ 

U.y4 / l\6l ) 


7 
z- / 


U.oUO 


in /Ifi 
1U.40 


1 Kfi/1 A 
1.0044 


1 good 
l.yzUU 




1.01 / / [Zi ) 


u.yoo / ^0 / ) 


7 
z- / 


U.oll 


11 /Ifi 
11.40 


1 fil Kf\ 

l.olOU 


z.Uzoo 




i fi7io/9/i\ 

1.0 / 1Z^Z4 ) 


u.youu^4i ) 


9 7 
Z- / 


n a 1 fi 

U.410 


1 9 A fi 
1Z.40 


1.0040 


9 1 Qfifl 
Z. 100U 


f 1 1 ft"* 


1 71/1/1/fiC\ 
1. / 144^0o ) 


u.oyoo^iou ) 


0-0 


n nnfi 
u.uuo 


10 /|(? 
10.40 


i 71 oa 
1. / 100 


9 9/191 
Z.Z4Z1 




i 77n / Q9\ 
1. / / Oi{6Z ) 


u.yio4(^oo ) 


9 £ 
z-o 


n con 
u.oyu 


1 /I A fi 
14.40 


1 7fi9f) 
1. / ozU 


9 Q/1 79 
Z.04 / Z 




1 QQfiO/QQ 1 * 

l.ooUz^oo j 


u.yiuy^oo ) 


9 

z-o 


n oqo 
u.uoy 


1 /I fi 
10.40 


1 qi nn 
l.olUU 


9 A K1 Q 
Z.4OI0 


(1,1a) 


l.o / / 0(^40 j 


u.oy4o(^ / 0) 


9 

Z-O 


n fico 
U.Ooy 


1 fi /I fi 
10.40 


l.oO / / 


z.OOOo 




i GQfifi/'/iG s \ 
l.yoUoi^y j 


U.00 / v{ IV) 


9 

z-0 


U.04o 


1 7 /I fi 
1 / .40 


1 OflKO 

i.yuoz 


9 fiKOp; 
z.ooyo 


1 ^ 


i QCfin/KfiA 
i.yoOU^OO J 


u.oooo^yu ) 


9 

z-0 


n fi/io 
u.04y 


10 /is 
lo.40 


1 QK9K 

i.yozo 


9 7fi9S 
Z. / OZo 


/"I 1 

(,l,loJ 


o nQ7Q/fi9 s \ 

z.Uo / o^oz j 


P) Q771 r oo^ 

u.o / / i^yy ) 


9 £ 

z-0 


U.1UU 


1 O A fi 

iy.40 


1 OOOfi 

i.yyyo 


z.oOOo 




i /] C71 /97\ 
1.40 / l^Z / ) 


n oo/i /i ^ 70^ 
U.yz44(^ / Zj 


Q fi 

o-O 


U.loU 


o.zO 


1 A 77Q 
1.4/ to 


1 QOQO 

l.oyoy 




i cn97/Q9\ 
l.OUz / [6Z ) 


u.yzyo^oy j 


O-O 


n 1 01 
U.lol 


Q O/l 


1 K991 

l.Ozzl 


1 ARKR 
1.4000 


(2,3J 


1.00lO^14j 


U.y440(^Z4 J 


9 /I 
Z-4 


U.Uoo 


O OQ 

y.yo 


l.OoUU 


1 KK7Q 

1.00 / 




1 fil CO/1 Q\ 

l.OlOz^lo ) 


fi QQQQ/ Q1 \ 

U.yooo^ol J 


9 7 
z- / 


n QQO 

u.ooy 


1 n oq 
iu. yo 


I.000O 


l.Oo / 




1 £71 Q/1 Q\ 

1.0 / lO^lo ) 


u.yoo / ^01 ) 


9 /I 
Z-4 


n /loo 
u.4yy 


1 1 

11. yo 


1 fiQ07 

i.oyu / 


1 7KQS 

1. / oyo 


(2,0) 


1 79/1G/99 s \ 

1. / Z4y^zz ) 


fl QQft^f QQ^ 

u.youo^oo j 


9 /I 
z-4 


n ofifi 
U.zOO 


1 9 OQ 

iz.yo 


1 7/1 Q1 
1. / 4ol 


i.oozo 




1 77Qn/97 s \ 

l.M oU^i / ) 


fl QO/I 1 / ^17^ 

u.yz4i ^4 / ) 


9 7 
z- ( 


u.4oy 


1 Q QQ 

10. yo 


1 7Q/1 1 

1. / y4l 


1 QfiKK 

i.yooo 


( 9 s\ 
(2,8) 


l.ooUO^zy ) 


fl Q1 fiO/ CfA 

u.yioz^ou j 


9 A 
Z-4 


U.U44 


1 /I OQ 

14. yo 


1 Q/1 A 1 
l.o441 


9 flfiQQ 

z.UOoo 


f9 c\\ 


l.o / oZ[oo ) 


fl QQ/I K/ KQ^ 

U.oy40(^0o J 


9 fi 
z-0 


9 9/1 1 
Z.z41 


1 K OQ 

10. yo 


1 fiOQQ 

l.oyoo 


9 1 7flfi 

z.l / Uo 


(2,10) 


f.9335(40) 


0.8985(68) 


2-6 


0.989 


16.93 


1.9420 


2.2732 


(2,H) 


1.9882(45) 


0.8961(74) 


2-6 


0.776 


17.93 


1.9902 


2.3755 


(2,12) 


2.0295(51) 


0.8681(82) 


2-5 


1.912 


18.93 


2.0381 


2.4776 


(2,13) 


2.0933(56) 


0.8818(92) 


2-4 


0.042 


19.93 


2.0857 


2.5796 


(2,14) 


2.1470(62) 


0.8768(102) 


2-5 


0.015 


20.93 


2.1331 


2.6815 


(2,15) 


2.1910(70) 


0.8552(115) 


2-7 


0.874 


21.93 


2.1804 


2.7833 


(2,16) 


2.2315(76) 


0.8272(119) 


2-5 


0.884 


22.93 


2.2274 


2.8850 
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TABLE IV: A part of lattice QCD results of the 4Q potential V±q for the symmetric twisted 4Q configuration labeled by (d, h). 
The notations are the same as in Table I. 



(d, h) 


ViQ 


C 


Tmm _ y max 


X 2 /Ndf 


L c4Q 
mm 


T/th 

c4Q 


V th - 

2QQ 


(3 1 1 


1 fifi41 (20) 


9086(341 


2-8 


1.052 


12.17 


1.7093 


1.5889 


(3 21 


1.6980(21) 


0.9095(35) 


2-4 


0.227 


12.65 


1.7359 


1.6314 


(3,3) 


1.7444(25) 


0.9098(43) 


2-6 


0.649 


13.42 


1.7769 


1.6941 


(3 41 


1.7960(26) 


0.9085(46) 


2-5 


0.278 


14.39 


1.8275 


1.7700 


(3,5) 


1 8473(31 "1 


f) 9D1 7(531 


2-5 


0.747 


15.39 


1.8787 


1.8540 


(3,6) 
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TABLE V: Lattice QCD results for the 4Q potential V4Q for the asymmetric planar 4Q configuration as shown in Fig. 6 with 
various (di,d 2 ,d s ,d 4 ) for h=8. The notations are the same as in Table I. 
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TABLE VI: Lattice QCD results for the 4Q potential Vaq for the asymmetric twisted 4Q configuration as shown in Fig. 7 with 
various (di, d 2 , d 3 , d 4 ) for h=8. The notations are the same as in Table I. 
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